Dynamics of Magnetic Defects in Heavy Fermion LiV 2 04 
from Stretched Exponential 7 Li NMR Relaxation 



in 
o 
o 



o 



i 

CO 

c3 



I 

C 

o 
o 



> 
m 

m 
o 
in 
o 

c3 



c 

o 
o 



X 
S3 



D. C. Johnston, S.-H. BaekQ X. Zong, F. Borsa, J. Schmalian, and S. KondcQ 

Ames Laboratory and Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011 

(Dated: March 21, 2005) 

7 Li NMR measurements on LiV204 from 0.5 to 4.2 K are reported. A small concentration of 
magnetic defects within the structure drastically changes the 7 Li nuclear magnetization relaxation 
versus time from a pure exponential as in pure L1V2O4 to a stretched exponential, indicating glassy 
behavior of the magnetic defects. The stretched exponential function is described as arising from a 
distribution of 7 Li nuclear spin-lattice relaxation rates and we present a model for the distribution 
in terms of the dynamics of the magnetic defects. Our results explain the origin of recent puzzling 
7 Li NMR literature data on L1V2O4 and our model is likely applicable to other glassy systems. 



Heavy fermion (HF, heavy Fermi liquid) behaviors 
have been widely observed at low temperatures T in 
many metals containing crystallographically ordered ar- 
rays of /-electron atoms, which are quite well understood 
theoretically 0. In these metals, the current carriers act 
as if they have a (heavy) mass that is of order 10 2 -10 3 
times the free electron mass. Only a few d-electron com- 
pounds are known to show HF behaviors at low T, e.g. 
Yi-zSocMns with x rj 0.03 2], LiV 2 4 00, and most 
recently C^-zSr^RuC^ with x ~ 0.3-0.5 [5j. There is 
currently no theoretical consensus on the mechanism for 
formation of the heavy fermion mass in L1V2O4 0. 

An important measurement for establishing Fermi liq- 
uid behavior at low T is nuclear magnetic resonance 
(NMR). For high magnetic purity samples of L1V2O4, 
the 7 Li nuclear spin-lattice relaxation rate 1/Ti is pro- 
portional to T (the Korringa law for a Fermi liquid) from 
about 10 K down to about 1.5 K 0, 0, 0| - In contrast to 
these results, recent 7 Li NMR measurements of several 
samples down to 30 mK by Trinkl, Kaps et al. strongly 
conflict with a Fermi liquid interpretation |9Lll0|. In par- 
ticular, non-exponential (stretched exponential) recovery 
of the nuclear magnetization, non-Korringa behavior in 
1/Ti versus T, a peak in 1/Ti at ~ 0.6 K, and a strong 
field dependence of 1/Ti were found at low T. In view 
of the small number of known d-electron HF compounds 
and the importance of LiV204 within this small family, 
it is critical to determine if these results are intrinsic to 
the pure material, and if not, what they are due to. 

Here we present 7 Li NMR measurements on two sam- 
ples from 0.5 to 4.2 K that were carried out to study the 
influence of magnetic defects on the low-T HF properties 
of LiV20 4 . We confirmed Fermi liquid behavior down 
to 0.5 K in a high magnetic purity sample. We find 
that a small concentration (0.7 mol%) of magnetic de- 
fects within the spinel structure drastically changes the 
detected spin dynamics and leads to the above behav- 
iors described in Refs. i and [H which therefore explain 
their results as arising from a significant concentration 
of magnetic defects in their samples. On the other hand, 
understanding the physics of magnetic defects in LIV2O4 
is interesting and important in its own right and may fur- 



ther guide and constrain theoretical models for the pure 
material. A crucial aspect of LiV204 is the geometric 
frustration of V spins for antiferromagnetic ordering in 
the spinel structure. The geometric frustration is likely 
directly related to the suppression of antiferromagnetic 
order in pure LiV204 and the emergence of a heavy elec- 
tron state instead. A large number of low lying spin ex- 
citations emerges and the system becomes "almost unsta- 
ble" , i.e. very susceptible with respect to crystal defects 
or other perturbations that locally lift the frustration and 
cause a condensation of the low lying states [llj| . 

We develop a phenomenological description for the ob- 
served stretched exponential 7 Li nuclear relaxation in 
terms of a distribution of 1/Ti values and explain the 
physical meaning of the parameters. We further present 
a model for defect nucleated dynamical magnetic order 
in an almost unstable electronic system that can explain 
this distribution and that provides important insights 
about the behaviors of magnetic defects in L1V2O4. The 
stretched exponential function also describes the kinetics 
of diverse relaxation phenomena [12j, so our model will 
likely have many applications to other fields. 

The two LiV204 samples measured here were samples 
#12-1 and #3-3-a2 studied previously in Ref. Hii where 
their synthesis and characterization were described. The 
magnetic defect concentration ridefect in the two samples 
was previously estimated from magnetization measure- 
ments at low T EH- Sample #12-1 shows a clear but 
weak intrinsic broad maximum in y(T) at 16 K, char- 
acteristic of high magnetic purity |3j, Il3| . with only a 
tiny upturn in x(T) below 4 K corresponding to ridefect 
= 0.01 mol%. The second sample, #3-3-a2, has ridefect 
= 0.7 mol% which is sufficiently large that the low-T in- 
trinsic broad maximum in x(T) is completely masked by 
the magnetic defect Curie-like term[l3j. The 7 Li NMR 
measurements were performed with a Fourier transform 
(FT) TecMag pulse spectrometer using 4 He (1.5-4.2 K) 
and 3 He (0.5-1.5 K) cryostats. The 7 Li NMR lineshape 
and the full width at half maximum (FWHM) were ob- 
tained from the FT of half of the echo signal. The 7 Li 
1/Ti was determined by monitoring the recovery of the 
spin echo intensity following a saturating pulse sequence 
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FIG. 1: (a) 7 Li NMR lineshape at 0.5 and 4.2 K and full width 
at half maximum peak intensity versus temperature (FWHM, 
inset) for two samples of LiV2C>4. (b) Hole-burning spectra 
at 4.2 K for sample #3-3-a2. 



of 7r/2 pulses. The typical ir/2 pulse width was 2 fis. 
The measurements were carried out at a frequency of 
17.6 MHz (magnetic field H = 10.6 kOe) so that a di- 
rect comparison of our results with the corresponding Li 
NMR data at 17.3 MHz in Refs.lil and Hoi could be made. 

The resonance line for the pure sample #12-1 has a 
FWHM that is independent of T below 4.2 K whereas 
sample #3-3-a2 has a much broader line [Fig. 1(a)] that 
becomes increasingly broad with decreasing T [inset, Fig. 
1(a)], indicating an increasing importance of magnetic in- 
homogeneity in the latter sample with decreasing T. The 
Knight shift K for both samples is the same and indepen- 
dent of T between 0.5 and 4.2 K, with K = 0.141(15)% 
in agreement with previous data |jj above 1.5 K. 

Our main experimental results originate from measure- 
ments of the influence of magnetic defects on the 7 Li nu- 
clear spin dynamics of LiV2 04- Figure 2(a) shows repre- 
sentative semilog plots of the time t dependent recovery 
of the 7 Li nuclear magnetization M(t) after initial satu- 
ration for magnetically pure sample #12-1 at several T. 
The data at each T lie on a straight line (shown) with 
a well-defined 1/Ti determined from a fit of the data by 
1 - M(t)/M(oo) = A exp^-t/Tx), where the prefactor 
A is typically 0.9 to 1.1. The resulting 1/Ti is plotted 
vs T in Fig. 3. These data follow the Korringa law for 
a Fermi liquid (1/Ti <x T) with a weighted fit giving 
(TiT)- 1 = 2.46(6)s" 1 K- 1 . Our results thus further con- 
firm Fermi liquid behavior for pure LiV204 at low T. 

The M(t) for sample #3-3-a2 with ridefect = 0.7 mol% 
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FIG. 2: Recovery of the 7 Li nuclear magnetization M vs time 
t for LiV 2 4 samples #12-1 (a) and #3-3-a2 (b). The solid 
curves in (a) and (b) are fits to the data at each T by an 
exponential and a stretched exponential, respectively. 
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FIG. 3: 7 Li nuclear spin- lattice relaxation rates in LiV2 04 vs 
temperature T. The lowest data set of 1/Ti vs T is for high 
purity sample #12-1; the fitted straight line is the Korringa 
law for a Fermi liquid. The 1/T X * vs T data sets are for our 
sample #3-3-a2 (filled squares) and from Kaps et al. (Ref. IToT) 
(open squares). The vertical arrow at about 0.7 K indicates 
the approximate T of the maxima in 1/T*. The inset shows 
the exponent f3 in the stretched exponential in Eq. (1) vs T. 



is shown for representative temperatures in Fig. 2(b). 
The recovery is drastically different from that of the 
pure sample, exhibiting strongly non-exponential behav- 
ior. Another important feature is that the shape of the 
recovery curve changes with decreasing T, particularly 
strongly below 1 K. Following Ref. Hfl we fitted the data 
at each T by the stretched exponential function 
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where (3 is the stretching exponent with < < 1 and 
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1/T* is a characteristic relaxation rate. This function 
is nonanalytic for t — ► 0. We therefore employed a low- 
t cutoff at 15 ms to the fit. The resulting 1/T-*(T) is 
plotted as hlled squares in Fig. 3 and the fits are shown 
by the solid curves in Fig. 2(b). The variation of (3 with T 
is shown in the inset of Fig. 3. The 7 Li NMR 1/T* data 
obtained at 17.3 MHz by Kaps et al. 01 are plotted vs 
T as open squares in Fig. 3. From the totality of the data 
in Figs. 1-3, we conclude that the data of Kaps et al. are 
not intrinsic to pure LiV204 but rather are dominated 
by the influence of magnetic defects. 

The non-exponential recovery of M(t) in sample #3-3- 
a2 suggests that there is a distribution of 1/Ti values in 
this sample for different 7 Li nuclei. To check this hypoth- 
esis, we performed relaxation measurements in different 
regions of the NMR spectrum. These are shown in Fig. 
1(b) as "hole-burning" experiments. By using a long sat- 
urating 7r/2 pulse we can irradiate only the central part of 
the spectrum. The fact that the "hole" disappears grad- 
ually during relaxation without affecting the remaining 
part of the line shape indicates the absence of spectral 
diffusion, which means that 7 Li nuclei with different Lar- 
mor frequencies have no thermal contact over our time 
scale. Another experiment was done by monitoring the 
recovery of the part of the echo signal far from t = 0. The 
recovery was found to be nearly exponential and with a 
1/Ti corresponding to the long time tail of the stretched 
exponential in Fig. 2(b). These two experiments together 
demonstrate that there does exist a distribution of 1/Ti 
values for the 7 Li nuclei in the sample at each T on our 
time scale and that these nuclei or groups of nuclei relax 
independently. The strong decrease of (3 with decreasing 
T in the inset of Fig. 3 must therefore reflect a signifi- 
cant change in the distribution of 1 /T\ values with T as 
discussed next. 

Our experiments thus demand that we model the 
stretched exponential relaxation in Eq. (1) as the sum 
over the sample of a probability distribution P of 1 jT\ for 
the various 7 Li nuclei. Accordingly we write the stretched 
exponential function in Eq. (1) as 

e-^ T '^ = / P{s,(3)e- st ' T *ds , (2) 
Jo 

where s — T*/T\ is the ratio of a particular relax- 
ation rate 1/Ti within the sample to the fixed relax- 
ation rate 1/T* characteristic of P(s,(3), and of course 
J °° P(s, (3)ds = 1. For (3 = 1, P(s, 1) is the Dirac delta 
function at s — 1. For general j3, Eq. (2) shows that 
P(s, (3) is the inverse Laplace transform of the stretched 
exponential. The evolution of P(s) with j3 for several 
values of (3 is shown in Fig. 4. With decreasing (3, P(s) 
broadens and becomes highly asymmetric, and the peak 
in P(s) becomes finite and moves towards slower rates 
which is compensated by a long tail to faster rates. The 
value of s at which P(s) is maximum is plotted versus f3 
in the inset of Fig. 4; this value decreases with decreas- 
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FIG. 4: Probability density P(s) for several stretched expo- 
nential exponents j3. The inset shows the s values versus j3 at 
which P(s) is maximum. 

ing (3 and approaches zero exponentially for (3 < 0.5. The 
physical significance of 1/Ti" f° r ?S < 1 is that 1/Ti 
is about equally likely to be less than 1/T* as it is to be 
greater; it is neither the average of 1/Ti nor the inverse 
of the average ofTi. For s » 1, P(s, /3) ~ l/s 1+/3 , so the 
average s avo = (T 1 */T 1 ) avo is infinite. Thus the moments 
of the distribution depend on the cut-off at large relax- 
ation rates, but this cutoff is irrelevant for the physical 
interpretation of the long time relaxation. We see that 
the measured small values for (3 at low T in Fig. 3 for 
sample #3-3-a2 and that of Kaps et al. constitute strong 
evidence for a broad distribution of 7 Li 1/Ti values at 
these T. 

We now discuss possible physical origins of the 1/Ti 
probability distributions in Fig. 4 and then propose a 
model that may be applicable to L1V2O4 containing mag- 
netic defects. We note at the outset that our NMR mea- 
surements were carried out at H — 10.6 kOe for which 
the magnetization of the magnetic defects at < 1 K is sig- 
nificant fl3| . a situation for which very few calculations 
of either the average bulk or local electronic spin fluctu- 
ations are available. With this caveat, the contribution 
to the 1/Ti of a nucleus located at site r, by fluctua- 
tions of electron spins S(r',£) with correlation function 
q r , y ,(t) = (S a (r'fi)S a (r",t)}, is 

^-(r)oc J d 3 r / d 3 r"A r yA r y,q r ^ r „{uj n ). (3) 

Here A r y is the hyperfine interaction between nuclear 
and electron spins at sites r and r', respectively, and 
Qr'.r" (w„) is the Fourier transform of the correlation func- 
tion q r ',r"(t) at the nuclear Larmor frequency u> n = j n H. 
For example, if one had a unique 1/Ti, a unique A r r i 
with r' = r", and a correlation function q r ', r '{t) ~ e~ £t , 
one would have 1/Ti °c e/(^ 2 + yielding a peak in 
1 /Ti as e decreases through uj n with decreasing T as ob- 
served for 1/Ti* at ~ 0.7 K in Fig. 3. This simplified 
example suggests that a significant fraction of the mag- 
netic defects drastically slows down below ~ 1 K. 
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A distribution of 1 jT\ could result from a spatial vari- 
ation in the electron spin dynamics, i.e. q r ',r"(t), or vari- 
ations of the hyperfine interactions. In the latter case, a 
nonlocal hyperfine interaction A r r > oc |r — r'| 3 , caused 
by dipolar and/or RKKY interactions, can lead to a 
broad distribution in 1/Ti. Depending on whether a 
given nuclear spin is close to or far away from a local de- 
fect that dominates the spin response q r ,r'(t), very differ- 
ent 1/Ti values occur. Geometric considerations lead to 
P(s, 13) oc s~ 3 / 2 , i.e. a fixed value = \ for the stretched 
exponential which is in direct conflict with our data that 
show a strongly T-dependent (3. A physically more inter- 
esting case is when stretched exponential nuclear relax- 
ation is due to dynamical heterogeneity of the magnetic 
defect spin system, i.e. due to spatially varying q r ' t r"(t). 
For simplicity we consider the limit of a purely local hy- 
perfine interaction A r r i oc 5 r - r i . In the limit of strongly 
disordered spin systems, dynamical heterogeneity with 
anomalous long time dynamics was found in numerical 
simulations above the spin glass temperature The 
averaged autocorrelation function was shown to have the 
Ogielski form q(t) = t~ x exp[— (e*t)@], where the energy 
scale e* characterizes the averaged electron spin response. 
As e* becomes smaller than u n with decreasing T, 1/Tj* 
goes through a maximum, located at T ~ 1 K in our case, 
and q (t) immediately yields a stretched exponential re- 
laxation for the nuclear spins. However, these results 
were obtained in the strong disorder limit in contrast to 
dilute magnetic defects in L1V2O4. 

The sensitivity of the HF state with respect to per- 
turbations may be critical to understand the behavior of 
L1V2O4 where ridefect is rather small, suggesting the fol- 
lowing alternative model for the dynamics. If geometric 
frustration suppresses long range order in pure L1V2O4, it 
becomes natural that crystal defects can locally lift the 
frustration and cause a condensation of dynamic mag- 
netic order in a finite region of volume ~ £ 3 around the 
defect. Due to the proximity to an ordered state and 
long-range electronic spin coupling in metallic systems, 
we expect £ to be much larger than an interatomic spac- 
ing, in contrast to insulating frustrated systems with only 
nearest- neighbor interactions ^^|. This might help to 
explain our previous low-T magnetization measurements 
which indicated that the magnetic defects have large av- 
erage spins ~ 3/2 to 4 13]. Fluctuations in the local ten- 
dency towards order lead to a probability p (£) oc e~ c ^ 
for such a droplet |lfi| and we analyze the system using 
the ideas of Griffiths physics in disordered magnets ^3] ■ 

The lowest excitation energies, e, of a droplet depend 
on its size £. Depending on how e (£) varies with £, 
different long-time dynamics emerges [18|. If e oc 
the distribution function of the droplet energies becomes 

p(e) oc e ~(^) 3/ "\ If 1/Tx oc e, this yields P(s,/3) oc 
e~ s , leading for large times to a stretched exponen- 
tial relaxation with (3 = 1+ ^y 3 - Such a behavior occurs in 



magnets with Heisenberg symmetry [18J where one finds 
/3 = |, if the spin dynamics is classical. For lower T, 
where quantum dynamics of the spins sets in, one finds 
(3 = |^-| < i, if the dynamical exponent relating length 
and time scales obeys z < 2. For insulating antiferro- 
magnets, one usually has 2 = 1, giving (3 = 1/3. 

Probably more appropriate to LiV204 is the case of 
itinerant antiferromagnets where one expects z = 2, 
and an even more exotic situation occurs. In this case, 
e oc e - b «'\ i.e. large droplets become extremely slow 
leading to quantum Griffiths behavior P (s) oc s~ x at 
long times with nonuniversal exponent A = 1 — c/b. 
Now, the nuclear spin relaxes according to a power law 
1 — M(t)/M{<x>) oc t~^~ x \ It becomes very hard to 
distinguish at large times power law from stretched ex- 
ponential behavior with small (3 at low T. However, 
there are clear predictions of this scenario which include 
(given A > 0) singular non- Fermi- liquid type specific heat 
C/T oc T~ A , susceptibility x K T~ x and similar results 
for the field dependence of the magnetization , which 
can all be tested in future experiments. 
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